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INCREASING STABILITY FOR THE CONDUCTIVITY AND 
ATTENUATION COEFFICIENTS 

VICTOR ISAKOV, RU-YU LAI, AND JENN-NAN WANG 


Abstract. In this work we consider stability of recovery of the conductivity and attenuation 
coefficients of the stationary Maxwell and Schrodinger equations from a complete set of 
(Cauchy) boundary data. By using complex geometrical optics solutions we derive some 
bounds which can be viewed as an evidence of increasing stability in these inverse problems 
when frequency is growing. 


1. Introduction 

The main goal of this paper is to demonstrate increasing stability of recovery of the conduc¬ 
tivity coefficient from results of all possible electromagnetic boundary measurements when 
frequency of the stationary waves is growing. When the frequency is large the Maxwell system 
can not be reduced to a scalar conductivity equation, so we have to handle the full Maxwell 
system. We do it by reducing the first order system to a vectorial Schrodinger type equation 
containing conductivity coefficient in matrix potential coefficient as in m, ns] and use ideas 
in m, m to derive increasing stability for potential coefficient involving variable attenuation. 
Use of the full system seems to be crucial, since increasing stability for the conductivity in 
the scalar equation does not seems to be true. We first consider a simpler case of scalar 
equation and extend the first and third author’s result [9] on the increasing stability estimate 
in the case of a constant attenuation to a variable one. Observe, that the method used in 
[9] works only for the constant attenuation. Extending this method to the case with variable 
attenuation is not an obvious routine, so the result for the scalar equation is of independent 
interest. 

The problem of recovering the conductivity from boundary measurements has been well 
studied since the 1980s. The uniqueness issue was first settled by Sylvester and Uhlmann in 
m where they constructed complex geometrical optics (CGO) solutions and proved global 
uniqueness of potential in the Schrodinger equation. A logarithmic stability estimate was 
obtained by Alessandrini [lj. Later, the log-type stability was shown to be optimal by Man- 
dache m ■ The logarithmic stability makes it impossible to design reconstruction algorithms 
with high resolution in practice since small errors in measurements will result in exponen¬ 
tially large errors in the reconstruction of the target material parameters. Nonetheless, in 
some cases, it has been observed numerically that the stability increases if one increases the 
frequency in the equation. The study of the increasing stability phenomenon has attracted 
a lot of attention recently. There are several results m, 0 , m and [T2i which rigorously 
demonstrated the increasing stability behaviours in different settings. 

For the Schrodinger equation without attenuation, the first author in [7] derived some 
bounds in different ranges of frequency which can be viewed as an evidence of increasing 
stability behaviour. The idea in [7] is to use complex- and real-valued geometrical optics 
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solutions to control low and high frequency ranges, respectively. The proof was simplified 
in [8] by using only complex-valued geometrical optics solutions. Recently, the first and 
third author showed in [9] that the increasing stability also holds for the the same problem 
with constant attenuation. The proofs in |9j use complex and bounded geometrical optics 
solutions. 

Continuing the research in [9] we show in this work even in the case of variable attenuation, 
the stability is improving when we increase the frequency. More precisely, the stability 
estimate will consist of two parts, one is the logarithmic part and the other one is the Lipschitz 
part (see Theorem EH]). In the high frequency regime, the logarithmic part becomes weaker 
and the stability behaves like a Lipschitz continuity. Moreover, the constant of the Lipschitz 
part grows only polynomially in terms of the frequency. We would like to point out that in 
m an increasing stability phenomenon was proved for the acoustic equation in which the 
constant associated with the Lipschitz part grows exponentially in the frequency. In order 
to obtain polynomially growing constants as in m , ei and [9], it seems like the smallness 
assumption on the attenuation is needed. 

The main result of the paper is a proof of the increasing stability for the Maxwell equations 
in a conductive medium. The first global uniqueness result for all three electromagnetic 
parameters of an isotropic medium was proved in [13] • We refer the reader to the latest results 
in [2] where the uniqueness and stability (log-type) were established using local data. Since 
our aim here is to demonstrate the increasing stability for the conductivity, we will assume 
that both the electric permittivity and the magnetic permeability are known constants. As in 
the previous results on the increasing stability, our proof here relies on CGO solutions. Here 
we will use CGO solutions constructed in [3]. Theorem 12.21 shows that the logarithmic part 
becomes weaker and the Lipschitz part becomes dominant when the frequency is increasing. 
In order to derive a polynomially frequency-dependent constant in the Lipschitz part, we 
impose the smallness assumption on the conductivity. Of course, it is an interesting question 
whether one can remove the smallness assumption and still obtain a polynomially frequency- 
dependent constant in the Lipschitz part. 

The paper is structured as follows. Section 2 introduces the mathematical settings for the 
Schrodinger and the Maxwell equations and presents our main results. The proof of stability 
estimates for the Schrodinger equation is carried out in section 3. In section 4 we derive the 
stability estimate for the Maxwell equations. 

2. Preliminaries and main results 

Let Q be a sub domain of the unit ball B in M 3 , with Lipschitz boundary dtt. By iL s (fi) we 
denote Sobolev spaces of functions on H with the standard norm || • ||( S )(H). These notations 
also hold for surfaces instead of H and for vector and matrix valued functions. 

2.1. Schrodinger equation with attenuation. We consider the Schrodinger equation with 
attenuation 

—Au — oj 2 u + qu = 0 in H, q = icua + c, (2.1) 

where <r, c are real-valued bounded measurable functions in fi. Since the boundary value 
problem for (12.11) does not necessarily has a unique solution, for the study of the inverse 
problem we consider the Cauchy data set defined by 

C{q) = ^(u,d u u)\dn € H l / 2 {d£l) x H^ 1 ^ 2 (dQ) : u is a IL 1 (H)-solution to (12.11) | . 
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Hereafter, v is the unit outer normal vector to d£\. Let qj = iujctj + Cj. To measure the 
distance between two Cauchy data sets, we define 


dist(C(<jfi), C{q 2 )) = sup sup inf 

j y k (fj,gj)dzCj (fk,9k)&Ck 


II (fj,9j) - (fk,9k) ||(1/2,-1/2) 

II (/.L 5.7 ) II (1/2,-1/2) (<9^) 


where ( fj,gj) = {uj,d u Uj),j = 1,2, are the Cauchy data for solutions Uj to the equation 
(12.11) with q = qj and 

1/2 


Denote 


II (/> 9 ) II ( 1/2 ,- 1 / 2 ) — (11/11(1/2) (<9f2) + || 5 ||(_i/2) 


e = dist(Ci,C 2 ), £ = — loge. 


With these notations, the bound which indicates the increasing stability phenomenon for 
the Schrodinger equation with attenuation can now be stated as follows. 


Theorem 2.1. Let s > 3/2 and qj = iuoj + Cj, j = 1,2. Suppose that supp(qi - ^) C 1! 
and Cj , (jj € H 2s {Ll) satisfy ||cj||( 2a )(fi) < M for some constant M > 0. Let u > 1,£ > 1. 

Then there exist constants m > 0 depending on s, M and C depending on s, D, M, such 
that if 

II a j II (2 S ) (^) < m, 

then 

II<?i — <? 2 1 |< C(u 2 e + + £) ^ s 3 ^ J ). (2-2) 

It is an immediate consequence of Theorem 12. II that if o\ = ci 2 = b in f l where b is a fixed 
constant, |6| < m, then one gets a bound similar to |9]. 

When | < s the logarithmic term is decreasing with respect to ui, so the bound (12.2p implies 
improving (better than logarithmic) stability of recovery of q from complete boundary data. 
We discuss it in detail. 

Obviously, 

+ w(w + S ) -(2»-3)/2 < w 2 e + (w + S) -(2s- 5 )/2 < w 2 £ + ^2 + £ 2) -(2 S -5)/4_ 

The function F(t) = et + (t + £ ,2 )-( 2s ~ 5 )/ 4 with the derivative F'(t) = e — (2s — 5)/4 (t + 
£■ 2 )—( 2 s—1)/4 - g i ncreas i n g Q n the interval [l,+oo) if a) (2s — 5)/4(l + £ 2 )~( 2s ~ 1 ')/ 4 < e and 
has minimum when (2s — 5)/4(t + £ 2 )~( 2s_1 )/ 4 = e if b) e < (2s — 5)/4(l + T 2 )^^ 2s_1 ^ 4 . 
In the exceptional case a) the minimal value of F on [l,+oo) is e + (1 + £ 2 )~( 2s_5 )/ 4 < 
e + (4/(2s — 5)e)( 2s_5 )/l 2s_1) . In the case b) the minimal value is 

e((4/(2s - 5)e)- 4 /( 2s - 1 ) - £ 2 ) + (4/(2 s - 5 ) e )( 2 ^-5)/(2s-i) < 


e((4/(2s - 5)e) -4 /( 2s-1 )) + (4/(2s - 5)e)( 2s - 5 )/( 2s - 1 ) = 

((4/(2s - 5)) _4 /( 2s_1 ) + (4/(2s - 5 ))( 2 *-5)/(2s-i) )e (2s-5)/(2s-i) i 

In the both cases, given error level e one can choose frequency u to guarantee at least Holder 
stability for q which is far better than logarithmic stability. Moreover, this analysis of the 
stability estimate (12.21) suggests that for given e there is an optimal choice of u> for the best 
reconstruction of q. 
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2.2. Maxwell equations. Let E and H denote the electric and magnetic vector fields in the 
medium with the electric permittivity e, p is the magnetic permeability //, and the electric 
conductivity a. The Maxwell equations at frequency u are 

f V x H + iu'yE = 0, , , 

\ V x E — iuipH = 0,7 = e + ia/uo. ^ 

We define the function space 

H (curl; fi) = {v € i/°(ft) : V x v € i/°(fi)}, 


with the norm 

IMI#(curi ; n) = IMI(o)(ft) + IIV x u||( 0 )(fi). 

For any v € //(curl; f2) the tangential trace of v, v x v, can be defined as an element of the 
space //~ 1 / 2 ((9fl). Namely, for any w* € i/ 1 (f/), 


{u x v, w) 

(see 0). We define the space 


/ (V x v) ■ w — 

Jo. 



X id), 


TH(dQ) = {w £ H l ^ 2 {d£t) : there exists d € //(curl; fi) , d 1 x v = w} 


with the norm 


IHIr.ff(dn) = inf{||w|U( CU r 1 ; n) : v x v = w}. 

The Cauchy data set corresponding to (12.311 is defined by 
C = {{y x E\ m ,v x H\gn) € ( TH(dQ )) 2 : 

(E,H) is (//(curl; fi)) 2 — solution of (|2.3p }. 


In this work we assume that e = p = 1. Suppose that Ci and C 2 are two Cauchy data 
of (12.31) corresponding to conductivities o\ and 07, respectively. We measure the distance 
between two Cauchy data C 1, C 2 of the Maxwell equations as follows: 


dist(Ci,C2) = sup sup 

j^k (F(j), G(j)) £ Cj , 


inf um,G(j)) - (F(k),G(km {TH{dn))2 

(F(k),G{k))eC k \\F(j)\\TH(dSl) 


where F{j) = v x E(j),G(j) = v x H(j) on d£l and E(j),H(j) is a H(curl;Q) solution to 
(12.311 with a = aj. Likewise, we denote 

e = dist(Ci, C 2 ), £ = — loge. 


Our main result is the following stability estimate for the Maxwell equations. 


Theorem 2.2. Let s > 3/2 be an integer. Suppose that supp(cri — 02) C 0. Assume that 
£ > c (depending on s,fi) and u > 1. 

Then there exist constants m and C depending on s,fl such that if 

IM(2a+2)(fi) <m ( 2 - 4 ) 

for j = 1,2, then 

||<7i — 02 II(_ s )(ft) < + £7 2 ) 3 / 2 e 1 / 2 + C(u + £)M 2s ~ 3 )/ 2 + C(uj + £)~ 1 . (2.5) 

As above, Theorem 2.2 implies increasing stability for recovery of the conductivity coeffi¬ 
cient. 
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3. Stability estimates for Schrodinger equation 

In this section we prove Theorem 12.11 We will use some ideas in [ 8 ] . C will stay for a 
generic constant depending only on s,Ll,M. First, following the arguments in [5], it is not 
hard to see that the following inequality holds. 


Proposition 3.1. Let u\ and U 2 be the solutions of \2.1\) corresponding to the coefficients 
(cti,ci) and ( 02 , 02 ), respectively. Then 

< II (ui, d v U\ ) 11(1/2,-1/2) || ( u 2i dvU2) II (1/2.-1/2) (<9^) e - (3-1) 


/ 

Jn 


( qi - q2)uiu 2 


We use CGO solutions for ( 12 . 111 . 


Proposition 3.2. Let £ € C 3 satisfy £ • £ = uj2 - Let s > 3/2. Then there exist constants Co 
and C\, depending on s and f l, such that if \C\ > Cb||<?||( 2 s)(f2), then there exists a solution 

u(x ) = e l< ^' x (l + ip(x)) 

of the equation \2.1\) with 

IMWfi) < ^lkll(2.)(fi). 

Proof. By extension theorems for Sobolev spaces there is an extension of q from Lipschitz Ll 
onto the unit ball B (denoted also q ) with |M|( 2 s)(-B) < C'(12)||q , ||(2 S )(^)- We can assume that 
q is compactly supported in B. 

We observe that since A)(e*^' x (l + if(x))) = (—A — 2 if • V + £ ■ £)(1 + f’(x)), if 

u = e z ^' x (l + if) is a solution of (12.11) if and only if the reminder if solves 

(—A — 2/£ ■ V)^ = — </(l + i/O in Q. (3-2) 

As known [15], there are constants Co, Ci, depending on s, such that if C'o 11 O' II ( 2 s) (-^) < |£|, 
there exists a solution if to (13.211 with 

M(2,)(B) < A-lkll(2,)(s). 

□ 


Let £ £ R 3 . We select e(j) € R 3 satisfying |e(l)| = |e(2)| = 1 and e(j) ■ £ = e(l) • e(2) = 0 
for j = 1, 2. Let R > 0. We choose 


1 


R 


2 \ 1/2 


2 , R 2 141A 1/2 


provided 


C(l) — — + * ( ~2~J e (!) + - 4 -J e ( 2 )> 

{(2) __*_ M *)%)_(^ 4 _]£)%), 


2 R 2 l£l 5 

w 2 + — > ^ 


It is clear that 

C(l) + C(2) = -£, c(i) ■ C(i) = w 2 , \((j)\ 2 = R 2 +u\ for j = 1 , 2 , 

which implies that \/2|£(j)| > R + uj. We observe that if 

R + uj > V2Co(muo + M), 


(3.3) 


(3.4) 
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then one has 

|C(j)l > C 0 (mu + M ) > C 0 \\iu(Tj + Cj||( 2s )(n). 

Note that when 

R > V 2 C 0 M and V2Com < 1, (3-5) 

the inequality (13.41) is clearly satisfied. So from Proposition 13.21 there exist CGO solutions 

Uj(x) = e jC(j) ' X (l + l/>j(x)), ||^j||(2 S )(0) < \\iuJCTj + Cy 11 (2s) (O). (3.6) 

To derive the stability, we would like to estimate the Fourier transform of q\ — q 2 - Before 
doing so, we need the following estimate on the Cauchy data of CGO solutions Uj, j = 1, 2 , 
constructed above. 

Lemma 3.3. If uj are the solutions LI. 01) . then 

|| (uj , dpUj )||( 1 / 2 ,— 1 / 2 )< Celi\/R 2 +uA (3.7) 

Proof. Recall that s > 3/2. Thus, by Sobolev embedding theorems 

lliMWn) < cUj\\ {2s) (n) < ^yWiuaj + Cj \\ {2s) (n) < < c 

and 

\ Uj (x)\ = | e <C(J>*(i + ^.( x ))| < e l»CC?)l(i + |^(x)|) < Ce^. 

Moreover, 

|Vu,(x)| < |iC0')e <C(7 ' ) -*(l + ^)l + k' <UVr Vr / | < Ce%(j)|(l + |W(*)|), 

hence 

||«j|| ( 1 ) (fi) < Ce^(l + |C(j)l)(l + ll^lld)^) < Ce^WR 2 + u 2 

and thus by Trace Theorems 

\\(uj, d u Uj)\\( 1/2 _ l/2) (dtt) < Ce75|C(j)| = Ce^\JR 2 Puf 2 . 

□ 


We denote / the Fourier transform of the function /. Let q = q\ — q 2 ■ Lemma [3.31 is used 
to bound the low frequency of q as the following lemma shows. 

Lemma 3.4. Let s > 3/2 and f = re with r > 0, |e| = 1. Assume that (13.51) is satisfied. 
There exists a constant C, such that 


|!(re)| < C(/>(R,u)Q(Q + Ce^(R 2 + w 2 )e, 
provided uj 2 + R 2 /2 > r 2 /A. Here 

Q(0 2 = J (^) _4S J(~i1 + f- x)~ 2s \q{ri)\ 2 diqdx 

and 


(3.8) 


<j)(R,u) = 


miw + M 


R + uj 
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Proof. Substituting CGO solutions (|3.6I) into (13.11) yields 

[ q{x)e~^' x (l + + ip 2 (x))d: 

Jn 


< 


||(^1; II (i/2,-1/2) (dty || (^2 j dpUf) II (1/2, —1/2) 


which leads to 

qe X dx 

q(x)e~^' x x(x)'h (x)dx 


I 

Jn 


< 


|| {ui, II (1/2,-1/2) (9fl) II (7^25 <9^7X2) II (1/2,—1/2) (9Ll)e, 

(3.9) 

where ’P(x) = i/q(x) + ^ 2 {x) + ifi(x)ip 2 (x), x € Cq°(Q), and y = 1 in supp(g). 

Since f fgdx = (27r) -3 f fgdg and fg = (27r)~ 3 / * g, one has 

= ( 2 »)-*/3(,) (e ^*)(-,)d, - (2 

Denote = e~ l ^' x x and (£) = (1 + |£| 2 ) 1//2 . By the Holder’s inequality, we obtain 

1 j 

< J \Q(v)\\(xi*^){-d)\dv 

< [ \q(v)\\ J Xs(~V ~ x)^(x)dx\d7j 

< I I \q(.v)\\xt(-V ~ x)\dri\if(x)\dx 


< Hx) 


—4s 


. _ \ 2 \ 1/2 
\q(.v)\M-'n ~ x)\drjJ dx j II^II( 2 s) 


<| / (x) 4s ( / (- ti + £-x) s \q{v)\(-V + f.-x) s \x{-g + f~x)\dg) dx ] ||^||( 2s ) 

\ 1/2 


1/2 


< (/( x ) AS (J(-V + t ~ x) 2s \q{g)\ 2 di^j dx^j 
Prom (13.61) and a priori assumptions, we have 


(s) II ^ II (2s) • 


11 t 11 ^ ^ mUJ + M 

11*11 P*)( fi)5C -RT7r' 

Therefore, the estimate (13.81) follows from (13.71) . (13.91) . (13.101) and (13.111) . 

The following lemma is an immediate consequence of Lemma 13.41 
Lemma 3.5. Let 

R* > V^CqM and \[2C§m < 1, 
where Co is the constant given in Proposition 14.41 


(3.10) 

(3.11) 
□ 


(3.12) 
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If 0 < r < lo + 1Z* , then 

|!(re)| < Cf{R*,u)Q(0 + C(lo 2 + R* 2 )e^ R * e; (3.13) 

if r > lo + R*, then 

\q(re)\ < C<j>(r,u)Q(£) + C(uo 2 + r 2 )e^ r e. (3-14) 

Proof. It is enough to take R = R* when 0 < r < lo + R* and take R = r when r > ui + R* 
in Lemma 13.41 □ 


With the help of Lemma 13.51 we can prove 

Lemma 3.6. Let the conditions of Lemma Iff.51 be satisfied. Then for any T > lo + R*, we 
have 

Il?llf_ 8) (n)<c'0(i2*, w ) 2 ||g||f_ 8) (n) 


+ C(lo 2 + T 2 )V^* + x{T)e 2V2T )e 2 + C{mu + m) 2 T" (2s " 3) , 
where x(T) < 1 and x(T) =0 ifT = lo + R*. 

Proof. Using polar coordinates, we obtain 

/*oo p 

||g||?_ s )(^) = / / |g(r , e)| 2 (l + r 2 )~ s r 2 dedr 

Jo J\e\=l 


pcu-\-R* p 

j'(re) 2 (1 + r 2 ) s r 2 dedr 

/ / 

Jo J |e|=l 

r T r 

+ [ [ 
Ju+R* J\e\ = l 

poo p 

g(re)| 2 (l + r 2 )~ s r 2 dedr 

+ / / II' 

JT J|e|=l 

: /1 + I2 + /:;■ 

(re)| 2 (l + r 2 )~ s r 2 dedr 


(3.15) 


(3.16) 


We begin with bounding I 3 . Since suppg C U by the Holder’s inequality |g(£)| < 
C 1 l<zll(o)( fi ) and so 

poo p 

I 3 <C\\q\\f 0) (n) / (1 +r 2 )~ s r 2 dedr 

Jt J |e|=l 

< C'||g|| 2 0) (H)T _ ^ 2s_3 ^ < C(uom + M) 2 T~( 2s ~ 3 \ (3.17) 

Before evaluating I\ and I 2 terms, we need the following estimate. Let A = { 77 ; | — i] + 
£ — x\ < |r?|/2}. By direct computation, we have 

J (0~ 2s J (x)~ 4s (j A {-V + €- x )~ 2B \Q('n)\ 2 d'nj dxdf, 

<cj <£>-/(x)" 4s (^j^{rf)~ 2s \q{r])\ 2 dr^j dxd£ 

<C'||g||?_ s) , 
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with C depends on s. On the other hand, by using the fact that A C {rj- ||£ — x\ < |ry| < 
2|£ — rc|} and (x)~ 2s (£)~ 2s < 2 s (x — £,)~ 2s , one has 

{—rj + £ - x)~ 2s \q(q)\ 2 dq] dxdf 
L / 

(—77 + £ — x)~ 2s \q(j])\ 2 drjdxd^ 


J ( 0 ~ 2s J {x)- is 

< J ar 2 ’ J (*)-“ 


4|£-®I<M<2|£-®| 


{x — Q 2s {-V + £,-x) 2s \q(q)\ 2 dqdfdx 


-2s | 


l£-*l<M< 2 |£-s| 


<C (x) 


< C 


-2s 


(v) 2s {-V + {,~ x ) 2s \q(r])\ 2 dr]d^dx 
(■x)~ 2s (-ri + £ - x)~ 2s d^dx] \q{Tq)\ 2 {rj)~ 2s dy 


l£-*]<M< 2 |£-x| 

-2s/ „ , a \-2s. 


<C\\q\\Uy 

Therefore, we can deduce that 


(0- 2S Q(0 2 ^<C\\q\\ 2 _ s y 


Next, by using (13.131) and (13.181) . we yield 

rui+R* 

h < / C((/)(R*,uj)Q( 0 + {u 2 + R* 2 )e V2Rt e) 2 (l + r 2 )~ s r 2 dr 

Jo 

<C(cj)(R\io) 2 \\q\\ 2 { _ s) {n) + (u; 2 + R* 2 ) 2 e V2R \ 2 ) J (1 + r 2 )~ s r 2 dr 


< Ccj>(R*,u;y\\q\\(_ s) (n) + C(u z + R* z ) z e 

since due to 2s > 3 we have 

puj+R' 


*2^2 a V2R* e 2 


rujy-rc no o 

/ (1 + r 2 )~ s r 2 dr < / (1 + r 2 )~ s r 2 dr = C. 

Jo Jo 

In the same way, the estimate (13.141) implies that 

h< C{cf(R*, u) 2 \\q\\ 2 { _ s] (n) + (u 2 + T 2 ) 2 e V2T e 2 ) f (1 + r 2 )~ s r 2 dr 


lu+R* 


<C<j>(ir,u)^q\\{_ a) (n)+C(u* + T*ye 


2\2 2\f2T 2 


e . 


(3.18) 


(3.19) 


(3.20) 


Combining (13.1611 - (13.201) and using that I 2 = 0 when T = u> + R*, we complete the proof. 

□ 

We are now ready to prove Theorem 12.11 

Proof of Theorem \2.1\ We will use the estimate (|3. 151) . Obviously, 

j. { 0 * n muj + M M 

<P(R ’ w ) = - <m + — ■ 

R* + oj R* 

Therefore, we can choose m < l/(y / 2C'o) and R* > \J2CqM so that 


c<kr*,u>y<-. 
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The choice of m depends only on s,D and of R* only depends on Thus, the first 

term on the right hand side of (13.151) can be absorbed by the left side, and we obtain 

||<?||(_ s) (ft) < C(oo 2 + r 2 ) 2 ^* + x {T)e 2V2T )e 2 + Coo 2 T~^ 2s - 3 \ (3.21) 

We consider the following two cases: 

(i) oo + R* < -£ and (ii) oo + 1Z* > —£. 

In the case (i) we choose T = We want to show that there exists C 2 such that 

{oo 2 + T 2 ) 2 e 2V2T e 2 < C 2 co 2 (co + £)" (2 *" 3) (3.22) 

and 

u 2 T -(2s-3 ) < c 2 w*{u + £)" (2 *~ 3) . (3.23) 

It is clear that (12.2|) is a consequence of 113.21(1 . 1(3.22(1 . and (13.231) . 

Due to our choice of T, (13.22(1 is satisfied if 

{00 + 8 fe 2 -^ 2 < C(co + £)~ 2s+3 , 

or 

{00 + £) a+2s ~ 3 e 2 ~ V2 < C, 

which follows from the elementary inequality (£) 4 + 2 s ~ 3 e 2_v/2 < (j 
Note that (13.231) is equivalent to 

c -l/{2 S -3)^ + £) < T = £ (3 24 ) 

Obviously, 

?>£ 

00 + £ < 00 + R* + £ < —, 

according to (i). Now we see that (13.241) holds whenever C;/ 1 ^ 25-3 ) <1/3. 

Next we consider case (ii). We choose T = 00 + R*. Then from (13.211) we have 

ll?lltf(- fl )( n ) < C (^ 4 + R* 4 )e 2y/2R *e 2 + Coo 2 (00 + R*)~( 2s ~ 3 \ 

Since R* < C and in the case (ii) 00 + R* > j + \£, the estimate (12.21) follows. 

The proof now is complete. 


4. Stability estimates for Maxwell equations 

In this section we will prove Theorem 12.21 We first discuss a reduction that transforms the 
Maxwell equations to the vectorial Schrodinger equation by adapting Section 1 in [4] to the 
special case of e = /r = 1. We then use arguments similar to the previous section to derive a 
bound of the difference of the conductivities. 
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4.1. Reduction to the Schrodinger equation. Let Q be a bounded domain in R 3 with 
smooth boundary. We consider the time-harmonic Maxwell equations (12.31) with e = fi = 1, 


[ V x H + iuj 7 E = 0, 
| V x E — iojH = 0, 


where 7 = 1 + ia/co, a € C 2 (fl) and a > 0 in Q. From (14.11) . one has the following 
compatibility conditions for E and H 


f V • ( 7 E) = 0, 

\ V ■ H = 0. 


(4.2) 


Let a = logy (the principal branch). Combining (14.11) and (14.2|) gives the following eight 
equations 

V • E + Va • E = 0, 

V x E — iojH = 0, 

' V • H = 0, 

V x H + iuj'yE = 0, 

which leads to the 8 x 8 system 



/ * 

0 * 

D- 



/ * 

0 

* Da- \ 



( 0 \ 


* 

0 * 

-Dx 


+ 

* 

ul 3 

* 0 



H 


* 

D- * 

0 


* 

0 

* 0 



0 


\ * 

Dx * 

0 

) 


\ * 

0 

* W 7/3 / 



V E / 


Here Ij is the (j x j)-identity matrix, 

D- = -i(di d 2 d 3 ), D = -i(di d 2 d 3 )* 

and 

/ 0 -d 3 d 2 \ 

Dx = —i I d 3 0 —d\ I . 

V ~ d 2 d\ 0 / 

* means that we obtain the same equations regardless of values of *. 
Define the operators 


P+(D) 


0 D- 
D Dx 


P-(D) 


0 D- \ 
D -Dx ) 


acting on 4-vectors. It is easy to check that 


P + (D)P.(D) = P.{D)P + {D) = -Ah 

and 

P+(D)* = P-(D), P-(D)* = P+(D). 
Also, we can see that the 8 x 8 operator 


is a self-adjoint, i.e., P* 


P = 


0 P.(D) \ 

PhD) 0 ) 


P. For C € C 3 , we set 


P( 0 = -* 


0 PhC) \ 
P+(0 0 J ' 
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Let X = (4>i, H t , $ 2 ) E t ) t be 8 -vectors, where $1 and $2 are additional scalar fields. The 
augmented Maxwell’s system is 

(P + V)X = 0 in SI, 

where 

( uj 0 0 Da- ^ 

0 01/3 Da 0 

0 0 017 0 

y 0 0 0 017/3 J 

Note that when <Lj = 0 the solution X corresponds to the solution of the original Maxwell’s 
system in f2. Now we want to rescale the augmented system. Let 

I4 . 0 

0 7 1/2 /4 


V = 


Y = 


then 


where 


(p + v)x=( 7 1 q Ia l 


X, 

(P + W)Y., 



( 0 

0 

0 

Da- 

\ 

1 

0 

0 

Da 

D x a 


2 

0 

0 

0 

0 



\0 

0 

0 

0 

/ 


with k = uiy 1 / 2 . Hence 

(P + V)X = 0 if and only if (P + W)Y = 0. 

The advantage of rescaling is that no first order terms appear in (14.3|i - (j4.5H . The proof of 
next lemma can be found in [4]. 

Lemma 4.1. 

(. P + W)(P- = -Ah + Q, 

(P - W l ){P + W) = -AI % + Q{ 1), 

(P + W*)(P - W) = -AI 8 + Q{2 ), 
where the matrix potentials are given by 


(4.3) 

(4.4) 

(4.5) 


( A a 


Q = 


1 


2 V 2 a — Aah 


\ 


( 1 


W = -2 


- K 2 I 8 ~ 

\ 


(Da ■ Da)h 


2 Dn- 


\ 2Dk 


2Dk 


2 Dk- \ 


(4.6) 


Aa 


- n z h - 


2 V 2 a — Aah ) 


and 


0 ( 2 ) = -- 


\ 


Aa 


- n 2 h - 


2 V 2 a — Aah ) 



\ 

2 Dkx 

—2 Dkx 

\(Da ■ Da)h 

( 

\ 


—2 Dux 

y 2 Dkx 

\(Da ■ Da)h 
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Here V 2 / = (d‘?j.f) is the Hessian of f and only non zero elements are shown in 8 x 8 matrices 
Q, Q(l), and Q(2). Moreover, W l denotes the transpose ofW and W* = W*. 

4.2. Construction of CGO solutions. We recall that C B. We extend a to a function 
in iL 2s+2 (M 3 ), still denoted by a, so that supper in B. Therefore, to 2 Is + Q, tv 2 Is + Q( 1) and 
uj 2 Is + Q( 2 ) are compactly supported in B. 

Let Y = (f 1 , ( u 1 ) t , f 2 , (u 2 ) t ) t €E H 2s (B). We recall the following construction of CGO 
solutions for Maxwell equations (14.11) in [3]. 

Proposition 4.2. [3l Proposition 9] Let s > 3/2 be an integer. Assume that a,b,( € C 3 and 
( ■ ( = lu 2 . Then there exists a constant Cq depending on s, such that if 


where 


ICI > Co ( ^ ||o ; 2 + Qj || (2s) (5) + ||w 2 /s + Qhs)(B) 
0 = 1,2 


q\ = —k 2 , q 2 = — -Aa — a 2 — -(Da ■ Da), 


then there exists a solution 


of (—A/g + Q)Z = 0 with 


Z = e i<: - x (A + ^) 




( C - a \ 

cob 


C ' b 

\ oja 


satisfying 


ll^ll(2 S )(fi)< ^\A\\\u 2 h + Qhs)(B). 

Moreover, Y = (P — W t )Z is a solution of (P + W)Y = 0 and 

Y = (0, H\ 0, 7 1 / 2 E t ) t 


where E, H solve (O- 

Proof. The proposition can be proved following [3] Proposition 9]. The only difference is that 
we use the estimate from H 2 f l to H 2s with — 1 < <5 < 0, that is, 

instead of [3 : , (22)] where the estimate is from L 2 +1 to L 2 , where G^ is the convolution 
operator with the fundamental solution of (—A — 2 if ■ V). This estimate (|4.7|) holds due 
to the easy fact (—A — 2 if ■ V) and (I — A) 2s commute. Due to the fact that u 2 Is + Q is 
compactly supported in B , one can replace the H 2 Y norm on the right hand side of (14.71) by 
H 2s (B). Recall that H 2s (B ) is a Banach algebra for s > 3/4. The rest of the proof follows 
[31 Proposition 9]. □ 


We also need CGO solutions for the rescaled system. These solutions can be constructed 
as in 131 Lemma 10, Proposition 11]. 
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Proposition 4.3. (3j Lemma 10, Proposition 11] Let a*,b*,( € C 3 and ( ■ ( = oj 2 . Then 
there exists a constant Co depending on s, such that if 

\(\>C 0 \\u 2 I 8 + Q(2)\\ {2 s)(B), 

then there is a solution 

Y* = e^ x {A* + tf*) 

of (P + W*)Y* = 0, where 

C • a * \ 

-C X a* 

C • b* 

c xb* 

and T* = PT* + iP( O'L* - WA* - W’tf* with A* = ^( 0 , b*\ 0 , a* 1 ) 1 , 

||**||( 2 S )(tt) < ^|^|||o ; 2 / 8 + Q(2)|| ( 2 s) (P) (4.8) 

and 

||P^*||( 2 S )(n) < Co(|A*| + ||'M (2 , ) (fl))|| W 2 / 8 + g(2)||( 2a) (5) (4.9) 

Moreover, 

ll*l (2s) (n) < Col^l (i + lkll(2 Y C | )(i?)a; ) (ll^s + Qm(2s)(B) + \\W\\( 2 s) (B)) . (4.10) 

4.3. Stability estimates. Suppose that e = p, = 1, ||<U ||( 2 s+ 2 )(-®) < m < 1 and supp(<Ji — 
a 2 ) C fh Hence we have supp( 7 i — 72 ) C fh We first prove an inequality that connects the 
unknowns and the boundary measurements. 

In the proofs C denote generic constants depending on s, H. 

Theorem 4.4. Assume that 07 and a 2 belong to H 2s+2 (Ll) and supp(<7i — < 72 ) C fi. There 
exists a constant C dependent on Ll,s such that, for any Z\ € H 2s (Ll) satisfying Y\ = (P — 
W{)Z\ = (0, H \, 0, r y\/ 2 E\) t with (Pi, Pi) solution to ^ with coefficient o\, and 

any H 2s (Q) solution Y 2 = (f 1 , ( u 1 ) t , / 2 , {u 2 ) 1 ) 1 of (P + W|)y 2 = 0, one has 

|((Q 1 -Q 2 )^ 1 ,l 2 )a|<C'dist(C 1 ,C 2 )||Ti||(i)(H)||r 2 || ( i)(H). (4.11) 

Note that Qj is defined in (14.61) corresponding to <jj for j = 1,2. The matrix-valued 

functions Wi, W 2 are defined similarly. 

Proof. Using supp(< 7 i — a 2 ) C H, from the first identity in the proof of (3) Proposition 7] one 
has 

((Qi ~ Q2)Z 1 ,Y 2 ) n = (Y 1 ,PY 2 ) n - (PY X ,Y 2 ) Q . (4.12) 

Moreover, using the estimate of the right hand of (14.121) derived in [3] (see the inequality 
before [3] Proposition 7]) and the support assumption of 07 — cr 2 , we can obtain that 

|(Yi,py 2 )n - (PYi,y 2 )n| 

< C(\\v x Ei — v x E 2 \\ TH ( d n) + \\u x H\ — u x ^ 2 ||rA/(av 2 ))II^II( 1 )(^) 

< C dist(Ci,C 2 )||y x ^i||TH(an)||> 2 ||(i)(n) 

< C'dist(C 1 ,C 2 )||yi|| ( 1 ) (H)||y 2 ||(i)(H), (4.13) 

where {y x E 2 , v x H- 2 ) € C 2 . The required estimate follows from (|4.12l) and (14.131) . □ 
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Let £ G M 3 . We select e(j) G M 3 satisfying |e(l)| = |e(2)| = 1 and e(j) ■ £ = e(l) ■ e(2) = 0 
for j = 1,2. Let R > 0. We choose 


2 7 


C(2) = -£-i 




where we assume 


Then 


2 R £ 2 

u 2 + — > — . 
2 “ 4 


C(l) - C(2) = -£, C U) • C(j) = w 2 , |C(j)| 2 =R 2 + w 2 , for j = 1, 2. 

To construct CGO solutions of Proposition 14.21 with Q = Qi corresponding to < 7 i, we 
choose a(l) = 0 and 

•»>~S*S 

By direct calculation, we can see that 

ll ^ 2 + Qj\\(2s)(R) + \\^ 2 h + Ql||( 2 -)(B) — C < 2 ||o'l||( 2 s+ 2 )(-S) w + C 3 , 
where C 2 , C 3 are independent of w. If we take 

R + uj > 2 1 / 2 (C 0 C 2 \\a 1 \\ { 2 s+ 2 ) (B)u + C0C3), 

then 

IC(l)l > 2 1//2 (i? + w) > Co (||uT + qj\\(2s)(R) + ||cD 2 i8 + (Jl|| ( 2 s) (-B)) • 

By Proposition 14.21 there exist solutions 

Z 1 = e i ^ 1) - X (A{1) + ^ 1 ) 

of (—A/g + Qi)Z\ = 0 with 


(4.14) 


/ 


A(l) = 


IC(1)I 


0 

ub( 1 ) 


\ 


C(l) ■ 6 ( 1 ) 

0 


and 


Co 


Likewise, if 


then one has 


|^l||(2 a )(fi)<|^Y)[l^(l)llK/8 + Q(l)||(2 a )(5). 

R + LO > 2 1 C(C'oC' 2 ||cr 2 ||( 2 s + 2 )(^)w + C 0 C 3 ), 


(4.15) 


|C(2)| >Co||( £ v 2 / 8 + Q(2)2|| (2 s)( j B), 

where Q(2) 2 is the matrix-valued function Q{2) corresponding to cr = er 2 . By choosing a* = 0 
and 

+ f 

in Proposition 14.31 with Q( 2) corresponding to < 72 , there exist solutions 

Y 2 = (A* (2) + %) 
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of (P + W%)Y 2 = 0 with 


A*( 2) = 


IC(2)| 


( 0 \ 

0 

C(2) • 6(2) 

V C(2) X 6(2) 


and + *P(C(2)) v h 2 * - W 2 A 2 * - W 2 \ h 2 *, where .4 2 * = |C(2)| 1 (0, 6(2), 0, 0) 

and if 2 * satishes (14.81) . (14.91) . Moreover, 

II ^2II (2s) (^2) < Co ^ ||cu 2 / 8 + Q(2) 2 || (2 s) (P) 


+ Co 


+ 


IC (2) | 2 ) 

W2\\(2s+2)( B )U 


IC(2)| IC(2)P 


I|w 2 || ( 2 S )( j b) 


(4.16) 


(see (14.101) ') . 

For Y 2 , we can obtain the following estimate. 


Lemma 4.5. Let Y 2 be the CGO solution of (P + 4F|)Y2 = 0 as above. Then for u > 1 there 
exists a constant C depending on s,Ll such that 

II V 2 || (1) {n) < C (1 + |C{2)|) e 2 - 1/a «. (4.17) 

Proof. Since V(e*C 2 )' x A* (2)) = i((2)e l ^ 2 ^' x A* (2) and |AL*(2)| < 2 1 / 2 , one gets 

11 ^(2)-*^*( 2 )|| (1) (^ < C (e 2 ~ 1/2 R \A* (2)\ + |C(2)|e 2_1/2jR |^(2)|) < C{ 1 + |C(2)|)e 2_1/2jR . 

By direct calculation, 

Wh + Q(2) 2 || (2 s) (P) < Cuj, \\W 2 \\ {2s) (B) < Cu, 

where we used the assumption 11 cr 2 1| ( 2 s + 2) (^2) <m<l. Thus, from (|4. 161) we have ||’l / 2 ||(i)(f2) < 
C 1 ^ 2 ) | which implies that 

||pC (2) -^*|| ( i)(ll) < || e ^ 2 )-*^|| (0) (fi) + ||<(2)e^ 2 >^|| (0) (tt) + || e <( 2 )-vvl/5|| ( o)(0) 


< C(1 + |C(2)|)^e 2 - 1/2 ^ < Cil + |C(2)|)e 2 1/2 ^. 


The proof is complete. 


□ 


Similarly, we can prove that 

Lemma 4.6. Let Z\ be the CGO solution of (—A/§ + Qi)Z\ = 0 and Y\ = (P — W^)Z\ as 
above. Then for uj > 1 there exists a constant C depending on s, LI such that 

ll*i||(i)( n ) < C(1 + |C(1)I) 2 e 2 ' 1/2H . (4.18) 

Proof. Substituting Z\ = (A{1) + ^i) gives 

Y 1 = (P- W{)Z i = (P - Wj)(e* (1 >- X (A(1) + Ti)). 

We evaluate 

< CICtOIII^WMfOUdjfS!) +C||ir; e i « 1 >M(l)|| (1) (S!) 

< C(1 + |C(i)l)V _,,lH 
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and 

IKP-^Ce^-^OUd)^) 

< C7|C(l)|||e^ (1) -*Wi||(i)(n) + C'||W^e^< 1 >*^ 1 || (1) (fi) + C F ||e i « 1 >‘ a! J P^ 1 || (1) (n) 

<c( i + |C(i)|)V 1/2R , 

which completes the proof. □ 

We now choose m < 1 satisfying 2 1 / 2 CoC 2 m < 1 and R > 2 1 ' /2 C'oC , 3, i.e., (14.141) and (14.151) 
hold. Combining (14.111) . (14.171) and (14.181) yields that 

\{{Qi ~ Q 2 )Zi,Y 2 )n\ < C dist(C 1 ,C 2 )||yi|| (1) (n)||F 2 || (1) (fi) 

< c( 1 + |C(2)|) 3 e 2l/2ii dist(C 1 ,C 2 ). (4.19) 


Lemma 4.7. Let s > 3/2 and f = re with r > 0, \e\ = 1. Under the assumptions of Theorem 
\4-4\ there exist a constant Cq depending on Ut only and a constant C depending on s,Ll such 
that 


1(5-1 - o- 2 ) (re) | < 


Cn + Cm„ fr , 2 ^ f ^ , C n m + Cm{ u;+ /?(£)) 


+ 


m2 -5(VV)(Q + 

Cnu; 1 /3(£) + Cm{u 2 + 


mf 


m 2 

Q(,m)+c 


s(ym) 


( a ; 2 + i ? 2 ) 7 / 2 2 1 / 2 


wm 2 


e 2 w dist(Ci,C 2 ), 
(4.20) 


where we denote 


d = a\ — a 2 , and 


m = R +[^ 2 + R 2 -^y\ ' , 


sum 2 = I (x)- 4 * /(->;+{- 


drjdx. 


Proof. Using Zi = e* < ’( 1 ^' :r (x4(l) + 4/i), Y 2 = e®£( 2 )' x (.A*(2) + and the definition of Qi, Q 2 , 
we have 

((Qi-Q 2 )Zi,Y 2 )n= [ e~^ x (AH 2 )+W 2 ) t (Q 1 -Q 2 )(A(l) + ^ 1 )dx 

JQ 

= f e~^ x A^2)\Q 1 - Q 2 )A(l)dx + [ e~ i ^ x ^f(Q 1 - Q 2 )A(l)dx 


+ / e-^ x {A*{2) + ^) t (Q 1 -Q 2 )T> 1 dx 

JVL 

= : h + I 2 + I 3 . 


(4.21) 



h 


(C(l)-6(l))(C(2)-b(2)) 

IC(1)IIC(2)| 



e 1 ^' x (k 2 — n\)dx 


+ 


2 ui 


C(2) - b( 2) 
IC(l)llC(2)| 


e 1 ^' x D(k 2 — k\) • 6(l)dx 


(4.22) 
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To bound I 2 , we recall that T?; = PT 2 * + zP(C(2))\I/ 2 * — W 2 A 2 * — W 2 'F 2 *, then 

h= [ e-^iP^YiQ, - Q 2 )A(l)dx + [ e-^ x (iP(C(2))^)\Q 1 - g 2 )^(l)dx 
Jn J n 

[ e^ x (-W 2 A 2 ^(Ch - Q 2 )A(l)dx + [ e-*- x (-W 2 * 2 *)\Qi ~ Q 2 )A(l)dx 
Jn Jn 


+ 


—: J\ + J 2 + J'A + Ja- 

We now bound Jk , k = 1,2,3,4. Applying (14.81) . (14.91) and using that m < l,u < |£(j)|, we 
can deduce that 


l*2*||(2 S )(fi) ga¬ 


rnet; 


IC(2)| 2 

By an argument similar to Lemma 13.41 we have 


, ||P^2*||(2 S )(^)<C' 


mu 

icM' 


(4.23) 


\Ji\ < c , Q(OI|P® 2 *ll (2 « ) ( n ) 
mul -Q((), 


< c 


IC(2)| 


where 

2(0 = 


IC(i)l 


<?(v 2 <x)(£) + 


w + C(l) -6(1) 
IC(1)I 


0(V<r)(O + 


cu 2 + w(C(l) • 6(1)) 
IC(1)I 


epm 


Using the bound ||P(£(2))\I/ 2 * ||( 2s )(fl) < Cmu |C(2)| 1 and following a similar argument i 


m 


Lemma 13.41 we can derive 


mu 


U 2 | < c^QK). 


For J4, it follows from the definition of W 2 that ||W 2 ||( 2s )(fi) < Cu, then applying (I4.23P and 
the proof of Lemma 13.41 we yield 

Finally, we would like to estimate J3. Note that 
W 2 A^ = 



f 

( 0 

0 

0 

Da 2 ■ 

1 

1 

K 2 h + ^ 

0 

0 

Da 2 

D x a 2 

IC(2)| 

0 

0 

0 

0 


V 

V 0 

0 

0 

0 


\\ 

l ° ^ 

6(2) 

0 

1/2 

^7 2 

1C(2) | 

O So 

-0 

// 

V 0 j 


V 0 / 


In view of the definition of Qj, direct calculations show that 
(Qi-Q 2 )A(1) 


(4.24) 


u 

( ° ^ 

2(V 2 (ai - a 2 ))6(l) - (A(«i - a 2 ))b( 1) 

iu(ai — a 2 ) 

( 0 \ 

ub( 1) 

2|C(1)I 

0 

V 0 ) 

IC(1)I 

C(l) - 6(1) 

V 0 / 


IC(1)| 


0 

- { u(Dai ■ Da\ — Pa 2 • Pa 2 )b(l) + 2(£(1) • b( 1))D(ki - k 2 ) 
2ub(l) ■ D(k 1 — k 2 ) 

0 


(4.25) 
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Now putting (|4.24j) . (|4.25l) together and using the fact that 6(2) • 6(1) = 0, we have 

(W2A2*Y(Qi — Q 2 )A(1) 

_ w 2 72 /2 6(2) t (V 2 (o:i - a 2 )) 6 (l) - 2w72 /2 (C(l) • 6(l))6(2) < H(/ci - re 2 ) 

IC(l)llC(2)| 

which implies, combined with to > 1, m < 1, |C(j)| > w, that 

|J 3 | < C n Qo(0 

where the constant Cq only depends on fl and 


Qo(0 = 


l 


IC(1)IIC(2)| 

Combining estimates for J\ to J 4 yields 


(atf(V 2 *)(0 + w(C(l) • 6(l))^(Va)(0 + (C(l) • b(l))G(am)- 


mix 


\h\<C—Q(ti) + CnQ 0 ((i)- 

To evaluate I 3 , from (|4.16l) . 

Il^ll( 2 .)(n) < ^1^(1)111^8 + Qi\\( 2 s)(B) < C- {1) + C(1) ' Kl) 


(4.26) 


|C(1)I 


ic(i)r 


it follows that 




(4.27) 


Note that ((j) ■ b(j) > 0 for j = 1,2. Combining (|4.19[> . (14.211) . (|4.22jl . (|4.26|) and (I4.27|) . we 
get 

(C(l)-6(1))(C(2)-6(2)) 


< 2lo 


IC(l)llC(2)| 
C(2) • 6(2) 
IC(1)IIC(2)| 


e~^ x {K 2 2 - K\)dx ) 

f e~ l ^' x D(n 2 — K\) ■ b(l)dx 

Jn 


mix 


+ Q(0 + CoQo(0 + F, (4.28) 

where F = (7(1 +|C(2)|) 3 e 2 1/2r dist(Ci,C 2 ). Substituting Kj = into the estimate (|4.28l) . 
we obtain that 

,(C(1) ' 6(1))(C(2) • 6(2)) 


ix- 


IC(l)llC(2)| 


/■ 


e + x (<72 — (j\)dx 


< 2ix 
which leads to 
|(di - 6-2)(01 < Cn 


[ e ^y 2 — 7 ^ 2 ) • b(l)da 

Jn 


C( 2) • 6(2) + C'^Q(O + ChQ 0 (O + +, 


1 


IC(1)HC(2)| IC (2) | 

,-l^ m lC(l)l 


+ C n 


(C(l) - 6(1)) 

IC(1)IIC(2)| 


(^(Va)(0 + w-W)(0) + c 


2o(0 + w 


(C(l) - 6(1))(C(2) - 6(2)) 
IC(1)IIC(2)| 


2(0 


w(C(l) - 6(1))(C(2) - 6(2)) 

IC(1)IIC(2)| 


< 91(0 + 92(0 + w " 


(C(l) ' 6(1))(C(2) • 6(2)) 


(C(l) - 6(1))(C(2) - 6(2)) 

F 


F 


(4.29) 
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with 

9i (0 = C n 
52 (0 = C 


1 


(C(l) ■ 6(1))(C(2) ■ 6(2)) 

m - 2 


G(V 2 am + Cn 


1 


(C(l) - 6(1)) 


G(Vam + Cn 


1 


«(C(1)-6(1)) 


gw® 


+ C 


(C(l) - 6(1))(C(2) - 6(2)) 
m(w 2 + w(C(l) • 6(1))) 


^(v a cr)(o + c g(Vff)(0 


(C(l) - 6(1))(C(2) - 6(2))' 


■ 0 (*)(fl- 


Since 


(C(l) ' 6(1))(C(2) • 6(2)) 

2 1/2 (C(1) • 6(1)) = 2 ~V 2 R + (<*> 2 + ^ ~ ' . (1 + IC(2)|) < 2(cn 2 + R 2 ) 1 / 2 , 


we have 

5i(0 < Ch^S(V 2 a)(0 + Ch^y0(V<7)(£) + 


(4.30) 


»(f) < cVL_S(vV)(f)+c ^y g^Kfl+(4.31) 


and 

where 


F < C(cj 2 + J R 2 ) 3 / 2 e 2 l/ 2 R dist(Ci,C 2 ), 


P^) = R+(2u; 2 + R 2 - l ^pj ' 


(4.32) 


The estimates (14.296 - (14.321) complete the proof. 


□ 


Recall that e = dist(Ci,C 2 ), £ = — loge. Before proving our main result, we need the 
following lemma. 

Lemma 4.8. Let R* > 2 1 ^ 2 CqCs with Cq and C 3 defined in Proposition \-j.3\ and (14.31) . 
respectively. There exist a constant Cq depending only on 0 and a constant C depending on 
s, 9 such that the following estimates hold: if 0 < r < oj + R*, then 


|(( 7 i - <r 2 )(re)| < 


Co + C' m ^/v 7 2 -uM , Cfi + Cm 


Mv 2 m) + 


R* + oj 


G(ym) 


(R* + oj ) 2 

+ ( ( f*\ , + Cm) G(m) + Coj~\oj 2 + R 2 *) 3/2 e 

\oj(R* + oj) J 


2*\3/2 0 2 1 t 2 R* £ . 


(4.33) 


if r > uj + R*, then 

|( 6 -i - o- 2 ) (re) | < 


Cn + Cm„ /r7 2~\ f/ -\ , Cn + Cm 


-g(v 2 dm + 


r + oj 


<?(vd)(0 


(r + oj) 2 

+ (+ Cm) G(a)(0 + Coj~\oj 2 + r 2 ) 3 / 2 e ‘^ 
\oj(r + uj) J 


(4.34) 


Proof. We choose R = R* when 0 < r < oj + R* and R = r when r >uj + R* in the estimate 

Ml . □ 

Now we are ready to show Theorem 12.21 
Proof of Theorem \2.2l Let a = a\ — er 2 . Using polar coordinates, we obtain 
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\a\\f_ s) (n)< / i*(oi 2 (i + i£i 2 r a de 

J\z\<U>+R* 

+ [ i^)i 2 a + iei 2 rx 

■J u+R* <\£\<T 
r oo 


r oo /• 

+ / / |cr(re)| 2 (l + r 2 )~ s r 2 dedr 

Jt J |e|=l 


l e l = 

= : 7 i + /2 + ^3 


(4.35) 


We begin with bounding J 3 . Since supper C 0 by the Holder’s inequality |<r(£)| < 
C||cr||(o)(71) < m < 1 and so 


h< 


poo p 

Cllall 2 (tt) / / (1 + r 2 )~ s r 2 dedr 

Jt J | e |=i 

< C7||cr||J 0) (^)T-( 2s “ 3 ) < CT-( 2s ~ 3 \ 


(4.36) 


Before evaluating ii and /2 terms, we need the following estimate. Let A = {t] : | — r) + 
£ — x| < |??|/2}. By direct computation, we have, for 0 < a < 2, 




—4 s 


(—V + £ — x) 2 s \V a a(r])\ 2 d'i] \ dxd^ 


< C J (0~' s J(x)~* s ) dxdt 

7Ct ~ 112 


< C'||V“cj|| 2 _ s) . 


(4.37) 


On the other hand, by using the fact that A C { 77 ; ||£ — x\ < \rj\ < 2|£ — x\} and 


(x) 2s (£) 2s <2 s (x — £) 2s , one has 

J (0~ 2s J (x)~ 4t 

< J®~ 2S JJ 2 _ 


f (—77 + £ — x) 2 s \V-a(ri)\ 2 drA dxd£ 


< C 

< C 


/«-//* 

I^iL 


3 l£-a|<H<2|£-a:| 


(—?7 + £ —x) 2 s \V a a(r])\ 2 dr]dxd £ 


(x — £) 2s (—77 + £ — x) 2 s \V a d(r])\ 2 dr]d^dx 

slf-z|<tol<2|£-x| 

(y)~ 2s (~V + £ — ic)^ 2s |V“u(77)| 2 (i77(i£(ix 

3l£-s|<H<2|£-a:| 

C f (/ J(x)- 2s (- V + C-x)- 2 s dCdx^ |V^(77)| 2 <77>- 2 s dr7 
< C||V Q d|| 2 _ s) . (4.38) 

We obtain from (14.37[> and (|4.38p that 


< 


J (£)- 2 ia(V“u)(£)| 2 d£ < C||V“d|| 2 _ s) . 


(4.39) 






22 


ISAKOV, LAI, AND WANG 


Next, by using (14.331) . (14.370 . (14.381) and (14.391) . we yield 

h< 


(C„ + c m f |e(v2 - )(s)|2 + M^, g(w)(f) p 


Cn 


|£|<uj+,r* ^ (R* + w) A {R* + u ) 2 

+ Cmj \G{d){f)\ 2 + Cuj~ 2 (uj 2 + R* 2 fe 2V2R \ 2 ^ {0~ 2s df 
|<r || 2 + Coo~ 2 (u 2 + R* 2 ) 3 e 2V2ir e 2 f 

‘ J|f 

+ Cw" 2 (w 2 + R * 2 ) 3 e 2y/ ^ R * e 2 . 


uj(R* + w) 


< Cm 2 ||d||^_ s) + 


C 


( 0 ~ 2 s d^ 


< C'm 2 ||(7||/_ Q x + 


{R* +lo ) 2 
Cm? 

'<-») ‘ r (i?* +w ) 2 
Similarly, we deduce that 

(Cn+Cl £t\g^am\ 2 + (c " +c " )2 |g(v^)(oi 2 


|£|<a;+.R* 


(4.40) 


I2 < 


L 


uj+R*<\£\<t'- (r + u) 


(r + oj ) 5 


+ f ^ . + Cm) \g(am\ 2 +u;- 2 (u; 2 + r 2 ) 3 e 2 ' /2r e 2 )(0~ 2s dt 

\oj(r + oj) J / 

< c™ 2 ll^ll?_ s) + (^pllSII^) + V + jf 


w+R*<|^|<T 


( 0 “ 2S ^ 


. 2 „.„ 2 , ^ + a ; - 2 ( c 4 ; 2 + r 2 ) 3 e 2 V 2 T e 2 _ 


(4.41) 


(4.42) 


< Cm ||d||f_ -, + 

" ll( sj (#* +w ) 2 

Therefore, we conclude from (14.351) . (14.361) . (14.401) and ((4.41 j> that 
Lemma 4.9. For any T > ui + R*, we have 

ll^„(«)< Cm’MU + J 0^ 2 

+ Coj~ 2 {uj 2 + T 2 ) 3 (e 2V2R * + X {T)e 2 V 2 T )e 2 + CT~? S ~ 3 \ 
where x(T) < 1 and x(T) = 0 if T = u + R*. 

We then choose m < 1 satisfying Cm 2 < 1/2. Thus, the first term on the right hand side 
of (14.421) can be absorbed by the left hand side. We obtain 

||*||(_ s) (n) < {r *° +uj) 2 + Cu~ 2 (u > 2 + T 2 ) 3 (e 2V2R * + X {T)e 2 V 2 T )e 2 + CT~? S ~ 3 \ 

Let 8 > 4 C 0 C 3 , we take 2y/2R* =8 so that R* still satisfies R* > 2 1 ^ 2 CqCs and we have 
||er||?_ 8 )(fi) < ° + Cuj- 2 {u : 2 + T 2 ) 3 {e~ l + X {T)e 2 V 2 T )e 2 + CT-? s ~ 3 \ (4.43) 

yO L>J J 

We consider the following two cases: 

(i) cj + R* < -8 and (ii) w + TV > —£. 

In the case (i) we choose T = §. By the same method as in the proof of Theorem 12.11 there 
exists C 4 such that 

u~ 2 (uj 2 + T 2 ) 3 e 2 V 2 T e 2 < C 4 (co + £)- (2 *" 3) (4.44) 

(4.45) 


and 


-V + T 2 fe 2 V 2 T e 2 < C 4 (w + £)" (2 *" 3) 
r -(2 S -3) < C 4 (w + 8 )~? S ~ 3 \ 
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Therefore, from (I4.43H - II4.45I) . one has 

11*11 (-.)(n) < + Cuj ~ 2(uj2 + f2)3e + c(w + V~ (2a ~ 3) - ( 4 - 46 ) 

Next we consider case (ii). We choose T = u + R* so that %(T) = 0. Then we have 
||a||J_ a) (fi) < + + Ccu~ 2 (uj 2 + (w + R*ffe + C{u + ip)- (2 *- 3) . 

Since R* = 2- 3 / 2 £, we have the same estimate d4.46D . We complete the proof of Theorem 

m 


5. Conclusion 

The next analytic task is to get rid of smallness assumptions of Theorems 2.1, 2.2 on 
attenuation and conductivity. It possibly will need a more careful reduction to a vectorial 
Schrodinger equation using some results in m It looks realistic to extend our increasing 
stability results for the conductivity coefficient onto the case of partial data [2]. To do so one 
can try to adjust the methods (in particular a new construction of CGO solutions) of the 
recent paper ma¬ 
lt is very important to verify if our theoretical prediction will result in better numerical 
reconstruction of the important conductivity coefficient for higher frequencies. After success¬ 
ful numerical testing one can expect a substantial improvements in the electrical impedance 
tomography, which so far was cheap, non destructive imaging technique with unfortunately 
pretty low resolution. 
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